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Abstract. We present a Rayleigh-Schrodinger-Goldstone perturbation formalism for 
many fermion systems. Based on this formalism, variational perturbation scheme 
which goes beyond the Gaussian approximation is developed. In order to go beyond 
the Gaussian approximation, we identify a parent Hamiltonian which has an effective 
Gaussian vacuum as a variational solution and carry out further perturbation with 
respect to the renormalized interaction using Goldstone's expansion. Perturbation 
rules for the ground state wavefunctional and energy are found. Useful commuting 
relations between operators and the Gaussian wavefunctional are also found, which 
could reduce the calculational efforts substantially. As examples, we calculate the first 
order correction to the Gaussian wavefunctional and the second order correction to the 
ground state of an electron gas system with the Yukawa-type interaction. 
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1. Introduction 

Field theories can be constructed from three kinds of pictures which are called 
Heisenberg, Interaction and Schrodinger pictures [1]. Among them, the Heisenberg 
picture is known to provide a convenient basis for the study of dynamics of operators 
and systematic perturbative improvement on physical quantities. On the other hand, 
the Schrodinger picture approach focuses on the dynamics of wavefunction, which allows 
detailed study on time evolution. However, in this approach, it is known that systematic 
improvement on obtained result is rather difficult. Instead, it has a powerful technique 
which is called a variational method which allows nonperturbative access to problems, 
so that it could be applied to strongly correlated systems where perturbative approaches 
break down. 

There have been successful applications of the Schrodinger picture to field theories 
[2-8] and later on to nonrelativistic many-particle systems [9-15]. While it has been 
proved useful, still the controlable range of trial wavefunctionals has been found to be 
narrow and further improvement from results of a trial wavefunctional demands much 
endeavor in field theories. Therefore, recent investigations on the picture have focused 
on overcoming its drawbacks towards a more managable theory. 

Gaussian trial wavefunctional apporach in bose fields has been widely applied and 
proven to be an efficient and powerful method in the Schrodinger picture [3,8]. Further 
improvements beyond the Gaussian approxiamtion have been investigated mostly in two 
directions. One is to try it with non-Gaussian wavefunctionals [16,17], and the other is 
to perform appropriate expansions based on Gaussian trial wavefunctional [18-24]. 

In fermi fields, a convenient prescription for operators was proposed by Floreanini 
and Jackiw [5]. It was shown that this approach gives successful results on Gaussian 
approximations of fermi fields as in bose fields. However, in contrast to the bosonic 
case [1], no successful Rayleigh-Schrodinger type perturbation formalism for fermi fields 
has been proposed so far, although several alternative schemes including functional 
integrals [25] and background field methods [8] have been reported using the Floreanini- 
Jackiw representation(FJR). Since the Rayleigh-Schrodinger perturbation formalism is 
the most familiar form of perturbation theories and especially suitable for the Schodinger 
picture representation in quantum mechanics, it is rather puzzling that it is not so in 
field theories. Indeed, the Rayleigh-Schrodinger perturbation scheme for bose fields was 
formulated early [1]. However, for fermi fields, Grassmann nature of functional variables 
has not allowed no such straighforward formalism as for bose fields. 

In this paper, we present a successful perturbation formalism for fermionic many- 
body systems combining the Rayleigh-Schrodinger perturbation with Goldstone's (RSG) 
expansion. The scheme will be used to formulate a variational perturbation scheme, 
which provides a systematic improvement beyond the Gaussian approximation. In order 
to make the presentation compact, we directly formulate a RSG perturbation scheme 
based on the variational Gaussian wavefunctionals. The formalism is trivially reduced 
to the conventional perturbation scheme by using a true Gaussian vacuum instead of an 
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effective Gaussian vacuum. Tliis point will be made clear once the formalism is given. 
It will be shown that the FJR provides surprisingly simple expressions for fermion 
wavefunctionals. As examples, we calculate the first order correction to the Gaussian 
wavefunctional and the second order energy of an electron gas with the Yukawa-type 
interaction. 

2. Floreanini-Jackiw representation 

The action of fermi field operators on Hilbcrt space can be represented by product or 
derivatives of Grassmann variables in the Grassmann function space. Then, Schrodinger 
equations arc transformed into functional differential equations [1]. Among many 
possible representations preserving the anticommuting relations of fermi field operators, 
the simplest one is generated from the fermion coherent states which are defined as 
eigenstates of annihilation operators [26]. In this representation, annihilation operators 
{aj} are described by Grassmann variables {u*}, and creation opeators {a|} , by 
Grassmann derivatives {^}- A representation without a complex conjugate notation 
* is also possible and has been used by Duncan et al. [4,9]. However, we find that the 
representation proposed by Floreanini and Jackiw is more convenient to formulate the 
present theory. 

In the FJR, the creation and annihilation operators arc given as follows [5]; 



where and u* are Grassmann variables. In the FJR for free fermi systems, vacuums are 
described by Gaussian wavefunctionals just as for bosonic systems. However, the norms 
between basis states for the FJR are not orthogonal since {u'u'*\uu* 
can be proven by the fact that four kinds of states relative to a quantum number i are 
possible. They are represented by 1, u*, Ui and u*Ui. Actually, only two kinds of physical 
states for i exist. One is the occupied state by a fermion and the other unoccupied. 
Thus, the FJR is a reducible representation and only a subspace of the total Grassmann 
functional space can be used to describe physically well-defined wavefunctionals. The 
subspace can be constructed only after a physical vacuum wavefunctional is chosen. The 
dual space wavefunctional ^ for the wavefunctional ^ is defined as 



where the norm is calculated by (^|^) = / L)[kk*]^^. It can be easily shown that the 
free fermion vacuum is described by Gaussian wavefunctional [5] as 



where G = and we find the dual wavefunctional is also Gaussian. Gaussian 

expectation values of some normal ordered operators are calculated as 




(1) 




(2) 



= DetG'^e^*^^"^ 



(3) 
(4) 




Rayleigh-Schrddinger-Goldstone variational perturbation theory for many fermion systems^ 

where / is an identity matrix and J, J* are source fields which are inserted into the 
Gaussian wavefunctional as e'^tj^t'^o+J'^t-'^tJi during the calculation. It allows that the 
Grassmann integrals are represented by Grassmann derivatives of source fields. 

For the free field Hamiltonian, Hq = hijalaj, the functional Schrodinger equation 
in the FJR is 

+ = ^»*- 

This equation contains second order derivatives and a quadratic term which are similar 
to the harmonic oscillator problem. Thus, we expect the Gaussian wavefunctional, 
Eq.(3), as a vacuum. In order for Eq.(3) to be an eigenfunctional with an eigenvalue 
£^0 — |Tr/i(/ + G) in Eq.(6), G should satisfy the following condition; 

(/ - G)h{I + G) = 0. (7) 

Eq.(7) has trivial and non-trivial solutions. Trivial solutions are G = ±/ and non- 
trivial ones, G = Here, we note that the number of particles N is given by 
N — ^Ttil + G). Therefore, G can be expressed as G = , where u is a 

chemical potential. Diagonalized, diagonal elements of G become 1(-1) below(above) /i. 
In the case of free Dirac fields, is zero and all the negative energy states are fully filled 
in the Dirac vacuum. 



3. Gaussian approximation and the state wavefunctionals 

When interactions between fermions exist, exact eigenfunctionals are different from 
the above Gaussian. Variational method or perturbation theory is applied in order to 
approximate true eigenfunctionals. Although the variational method depends largely on 
intuition in contrast to the systematic perturbative approach, it could provide excellent 
results if trial states are chosen carefully. 

In the Gaussian approximation, Gaussian wavefunctional with variational 
parameters is used as a trial functional. We now consider the case of interacting fermion 
system where the Hamiltonian has a general form, 

H = hijajaj + Vij^aja^akai. (8) 

Trial wavefunctional is chosen to be 

^g = e^*^-<"^ (9) 

where G is a variational parameter matrix. Then, the energy expectation value is readily 
calculated with the aid of Eq.(5). 

where Sj^ = {vkiji — Vikji)^ik and Q = {I + G){G + G)^^{I + G). Minimization 
of the energy under a fixed number of particles is achieved through the relation. 
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E — n(^Tril — N) = 0, where fx, the chemical potential, is introduced as a 
Lagrange's multiplier. Thus, the solution G should satisfy the following condition; 

(7 - G) (/i + E - (7 + G) = 0. (11) 

The equivalent condition for G is obtained by minimizing with respect to G, and is 
given by {I + G) {h + T, — jil) [I — G) = which provides the same G as in Eq.(ll). As 
in Eq.(7), the nontrivial solution G of Eq.(ll) is given by 

G = -^i±£^^. (12) 
^(/i + E-/.7)2 

With this solution, the ground state energy Eq is given as Eq = |Tr(/i + S)(7 + G). 
The parent Hamiltonian 77o which has the above Gaussian vacuum is easily found to be 
^0 — {h-\-T)ija\aj. In order to obtain excitations based on the Gaussian wavefunctional, 
we should have a unitary matrix U which diagonalizes G as UGW — 7, where 7 is a 
diagonal matrix with elements of +1(-1) below(above) the Fermi level Thus, G can 
be rewritten as follows; 

Gi G2 
G2 G4 

Therefore, we obtain an expression for U, 




The Gaussian wavefunctional is diagonalized by the new Grassmann variables Ui — UijUj 



and u* = u*U'ji. We use capital(small) letters above(below) the fermi see, so that 



We define new creation and annihilation operators as follows; 

Excited wavefunctionals and their duals are obtained by these operations to the Gaussian 

as 

^excited ^U^--- uXUa, ■ ■ ■ Ma„e^»^^>^ , (17) 
^excited ^Kr.--- K.'^Ar. ' ' ' ^^Aie^'^^"^' , (18) 

which has an excitation energy of e^j + • • • + ca^ — (^ai + ■ ■ ■ + ^a„) where ej is the 
ith matrix element of U{h + S)?7^ We note here that any physical wavefunctional in 
the FJR could be represented by multiplying only some combinations of u\ and Ua 
to the Gaussian, while in other representations, wavefunctionals should also contain 
Grassmann derivatives for describing hole states and, thus, arc quite complicated in 
general. An important aspect of the present result is that Eqs.(17) and (18) have same 
structures to those in bose systems [1] . 
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4. Rayleigh-Schrodinger-Goldstone Perturbation Formalism beyond the 
Gaussian approximation 

In order to go beyond the Gaussian approximation by a perturbative method based on 
the parent Hamiltonian Hq, Hamiltonian is rearranged as 

H = Ho + {H- Ho) 

= {h + S)ija|aj + (Vijkiala]akai - Sijafaj) 

^{h + T,)ijaldj + (Vijkidlalakai - ^ijdldj^ , (19) 

where h — UhW, E = UHW and Vijki — UuUjmVimpgUpk^li- Here, we note 
that the whole formalism simply reduces to the conventional perturbation if E is 
set to zero. In the previous section, we have already obtained whole spectrums 
of eigenvalues and eigenfunctionals of the parent Hamiltonian Hq. Therefore, the 
Raylcigh-Schrodinger perturbation procedure [1] can be readily adopted for improving 
the Gaussian approximation. However, in order to carry out the calculation, it is 
necessary to have a general rule for combination of Grassmann operations which is 
similar to Wick's theorem of the Green function approach. In the following, we show that 
another perturbative approach in the time independent formulation, namely Goldstone's 
expansion [27], becomes an extremely useful tool in the present approach. Since final 
results from Goldstone's expansion are same as from the Raylcigh-Schrodinger scheme 
order by order, we follow the Goldstone's approach to obtain the expressions for the 
Rayleigh- S chrodinger per t ur bat ion . 

Goldstone's expansion is as follows. U H = Hq + Hi, Hq\^q) = Eq\^q) and 
i/|*o) = {Eo + AE)\^o), then 

l*o) = E'(;^— ^^i) I$o)l, (20) 
/I \ " 

AE^J2'{^o\Hi(- -Hi) |$o)c, (21) 

where |^o) is normalized as ($o|^o) = 1- The prime in summation represents exclusion 
of terms with zero denominators and subscripts L and C represent linked and connected 
diagrams respectively. The connected has the same meaning as in other many-particle 
theories. However, the linked diagrams in the Goldstone's usage are such that even 
when some diagrams have unconnected parts, if unconnected parts all have external 
lines, then, the diagrams are treated linked. 

Let's expand the exact ground state wavefunctional ^ in terms of the renormalized 
interaction, H - Hq. It gives * = En*n, where *o = *G = e*^*^"*"^. Then, the nth 
order wavefunctional is given by 

*n = E' i-^r^iH - Hq)] e^^^<-^, (22) 
where / D[u*u]^g^ = / D[u*u]^g'^g is satisfied. 
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In order to represent explicitly the Grassmann wavefunctional, we use the following 
useful relations. When capital(small) letters are defined for states above(below) the 
fermi see, that is, A,B,--- > kp and a,b,--- < kp, we arrive at simple commuting 
relations as follows, 

<*G = *G- 



aa'^G = ^GV^«a. (23) 

These commuting relations will be shown to reduce calculational efforts substantially in 
the FJR. We further define i and i for notational simplicity as 

dl^G = ^Gi, 

ai^G = *G«, (24) 

where explicit forms of i and i in the right-hand side are defined through Eq.(23) 
according to their position in the fermi see. Then, we obtain 

*n = E'e''^^"'*"^ (tt^ {V^jkimi) - ^^J&))] ■ (25) 

V-fc/0 - -Ho J L 

The right-hand side above consists of combinations of Grassmann variables and 
Grassmann derivatives. We can easily see that if there is no Ui (u*) on the right-hand 
side of ^ (af^)) it should be zero. Contraction is defined as a Grassmann derivative 
result which is expressed by 

S = £-uje{kF - i) = 5ije{kF - i), (26) 

?j = i^upii - kp) = SijOii - kp). (27) 
Self-contractions of bare interaction, Vijki{ijkl) are canceled out by T^ij^ij) because 

= {Vhji - v,,ji){unWMtj) = 2{v,,ji - v,kji)6ik{e{kF - /) - e(/ - kpMj). 

Therefore, the interaction term can be simply expressed as 

%ki{ijkl) - tij{tj) = %u{f3kl)\ (28) 

where the prime means that contractions should be performed with elements outside 
{ijkl). With this prime notation, the wavefunctional is expressed as 



^ n 



6Ei 6E2 5En 

X [{{iljlkili)' {i2j2k2l2)' ■ ■ ■ {injnknln)')!] , (29) 

where 6E^ = Yla=^{^ia + ^ja ~ ~ ^l) means excitation energy. Each index in the 
second square bracket in Eq.(29) could be a Grassmann variable or a Grassmann 
derivative. Therefore, in order to be non-zero, all derivatives should be contracted 
out with corresponding variables and finally remaining elements should be Grassmann 
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variables only. They are summarized as follows: (i) The following steps arc repeated for 
each m = [0, 2n — l]. (ii) A possible set of m pairs of bar and non-bar indices is selected. 
Self-contractions are forbidden and isolated parts fully connected by contractions are 
also excluded, (iii) Each contraction has the value given by Eq.(26) or Eq.(27). If odd 
numbers of permutations are performed for contraction, (-1) is multiplied. Remaining 
indices without contraction are converted by i — > \/2u*Q{i — kF) and i \/2uiQ{kF — i)- 
(iv) All other possible m contractions are performed. Then, is calculated by 
multiplying the first square bracket term of Eq.(29) to each contraction term and finally 
by summing out all indices. 

The corresponding nth order energy £^„ {AE — J2n ^n) is calculated as 

_ !D[u*u\^G{H-H^)-^n 
JD[u*u]^g'^g 

E't> ~^2i2fc2^2 ~'^ininknln 

X {{iljlkih)' {i2j2k2l2y ■ ■ ■ {in3nknln)')c > (30) 

where subscript C means that we should perform fully connected contractions. 

As well as energy values, Gaussian expectation values of any operators like 
0{a\^ ttj, ■ ■ •) are also simply expressed in the FJR as follows. 



J D[u*u]e^^^^''>iO{i,j,- ■ •) 



= 0{Ij,---)\u',u=o, (31) 

where G should be the variational solution, Eq.(12), and the notational abbreviations, 
Eq.(16) and (24) are also used. The third line of Eq.(31) is equivalent to the full 
contractions allowing unconnected ones. 

As examples, we calculate the first order wavefunctional ^'i and the second order 
energy £2- If we use capital(small) letters above(below) the Fermi see in the correction 
parts, they are represented as 

^^ = _4eG..<«. E ~ ^«»^a«6, (32) 

A,B,a,b ^A + eB-ea- 6^ 

£i = 0, (33) 

g _ _ ^VabAB {VABab — VsAab) ^^^-^ 
^A + ^B — — h 

Here, we note that £i is always zero and has no terms hke u\Ua because self- 
contraction is forbidden. 

In the case of electron gas, the quantum numbers are usual wavevector k and spin 
a and the interaction has a form, 

Vijkl ^ ^fcicri,A;2Cr2,fc3(73,fc4(74 = o'^(l^4 ~ kl\)Ski-\-k2,k3+k4^iJia4^a2CT3 ■ (35) 
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For a homogeneous case, the Gaussian variational solution G of Eq.(12) is already 
diagonalized by ka representation. Therefore, we can use the Grassmann variables 
■u^^, Uka as basis states instead of u\^, Uka- The energy component of parent Hamiltonian 
Hq is calculated as ika = [h + Ti)ka,ka = ^k ~ J2qi^{(l)''^k+qo- where e° is the free electron 
energy band and rika is the Fermi distribution function at zero temperature. Thus, we 
finally have 

^ _2e«-.'=v-L-.v' Y: z ::-ul,_,^,ul_^,,Uka,Uk>., (36) 

kk'qaa' ^'='+<? + ^^-q - Cfc - Cfe' 
C '^{QynkTlk'il - nfc/+q)(l - rik-q) 

^2 - - 2^ - 



kk'q 



^k'+q + ^k-q — ^k — ^k' 



^ v{qMq')nk+q+q'{l - nk+g')nk{l - Uk+q) ^2^^ 

kk'q ^k+q + ^k+q' ~ ^k ~ ^k+q+q' 

The first order wavefunctional has a very simple form and higher order ones are 
also expected to have such simple forms in this variational FJR. We find that the 
above ground state energy is same as the zero temperature limit of the variational 
perturbation theory [25] previously obtained from a functional integral representation, 
because Gaussian approximations of both representations employed same variational 
basis. For a numerial energy value, we modeled an electron gas with a bare kinetic 
energy = and a Yukawa-type interaction v{q) — Ai^'-^e^ /V{q^ + {^Y) where me 
and e are the bare electron mass and the charge, and V and Oq are the system volume 
and the Bohr's radius respectively, a, 7, and A are parameters to modulate the mass, 
charge and interaction range. In Fig.l, we plot the second order contribution to the 
ground state energy as a function of r^, which is the average distance ratio between 
electrons defined hj V = ^ti {r sOq)'^ N . We find the second order contribution is much 
reduced in the present variationl procedure than the conventional one which usually gives 
exceedingly negative result. This indicates much better convergence for the variational 
perturbation theory as expected. 



5. Summary 

We have presented a Rayleigh-Schrodinger-Goldstone perturbation formalism using the 
Floreanini-Jackiw representation on fermion systems. With the aid of Goldstone's 
expansion, formal expressions for the ground state wavefunctional and the corresponding 
energy are obtained in terms of a renormalized interaction. It is shown that this 
formalism can be used for both standard perturbation and variational perturbation 
schemes. Useful commuting relations between creation and annihilation operators with 
the Gaussian wavefunctional have been found. As examples, we have calculated the first 
order wavefunctional and the sencond order ground state energy of electron gas with 
Yukawa- type interaction. 
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Figure captions 



Figure 1. The second order energy correction (Eq.(36)) per particle for the electron 
gas with the Yukawa-type interaction, (a = 1.1, 7 = 1.1 and A = 0.1) ; filled circles 
are from the present variational perturbation and open circles from the conventional 
perturbation. 
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